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SUMMARY 

A first-principles theory is developed to study the nonlinear spatial evo- 
lution of a near-resonance triad of instability waves in boundary- layer transi- 
tion. This triad consists of a plane wave at fundamental frequency and a pair 
of symmetrical , oblique waves at the subharmonic frequency. A low-frequency, 
high-Reynolds-number asymptotic scaling leads to a distinct critical layer 
where nonlinearity first becomes important; the development of the triad s 
waves is determined by the critical layer's nonlinear, viscous dynamics. The 
resulting theory is fully nonlinear in that all nonl inear ly generated oscilla- 
tory and nonosci 1 latory components are accounted for. 

The presence of the plane wave initially causes exponent ial-of-exponent ial 
growth of the oblique waves. However, the plane wave continues to follow the 
linear theory, even when the oblique waves' amplitude attains the same order of 
magnitude as that of the plane wave. A fully interactive stage then comes into 
effect when the oblique waves exceed a certain level compared to that ot tne 
plane wave. The oblique waves react back on the fundamental, slowing its 
growth rate. The oblique waves’ saturation results from their self-interaction 
- a mechanism that does not require the presence of the plane wave. The 
oblique waves' saturation level is independent of their initial level but 
decreases as the obliqueness angle increases. 

Part II of this study will present results for the composite solution and 
comparisons with observations. Novel features of the phenomena will be 
described, and new interpretations of the experimental data will be given 
therein. 


1. INTRODUCTION 
l.l. Background 

Experiments on boundary- layer transition and accompanying instability-wave 
growth have identified several distinct flow-development regimes. The first 
regime begins with the onset of a two-dimensional, Tollmien-Schlichting (TS) 
wave that propagates in the flow direction and is uniform in the spanwise 
direction. These waves harmlessly grow and decay in accordance with the two- 
dimensional, linear stability theory when their amplitudes are sufficiently 
small. However, a nonlinear regime characterized by a three-dimensional dis- 
turbance field appears when the amplitude becomes sufficiently large. These 



three-dimensional disturbances grow at much higher rates than predicted by 
linear theory and ultimately lead to the transition to the turbulent state. 

The occurrence and growth of three-dimensionality have been recognized as 
a prerequisite for the eventual transition into turbulence since the experi- 
ments of Klebanoff & Tidstrom (1959) and of Klebanoff, Tidstrom, & Sargent 
(1962), which demonstrated that the three-dimensional structures are periodic 
in the spanwise direction. These experiments, which incorporated relatively 
large two-dimensional input disturbances, documented the wrapping of the ini- 
tial TS wave into peak-valley pairs characterized by streamwise wavelengths 
equal to that of the TS wave and to a spanwise wavelength of the same order. 

But at low-to-intermediate levels of the two-dimensional input disturbances, 
the experiments of Knapp & Roache (1968); Kachanov, Kozlov, & Levchenko (1978); 
Kachanov & Levchenko (1984); Saric & Thomas (1983); Saric, Kozlov, & Levchenko 
(1984); and Corke & Mangano (1987, 1989) indicated the existence of another 
form of the spanwise periodic structure. In these experiments, a staggered 
X-shaped pattern occurred, the streamwise wavelength was twice that of the TS 
wave, and the spanwise wavenumber was almost half of that observed by Klebanoff 
et ai. The three-dimensional structure is thus the subharmonic of the input TS 
wave. Unlike the TS wave, the observed subharmonic oblique waves grow rapidly 
and eventually lead to the transition to the turbulent state. Since this struc- 
ture occurs at low-to-intermediate amplitudes of the TS wave, it is more likely 
to occur in natural transition than is the fundamental- type structure of 
Klebanoff et ai . These experimental observations established the subharmonic 
route to boundary- layer transition, which is the subject of the present work. 

The observed three-dimensionality in boundary- layer transition has promp- 
ted several theoretical and numerical attempts to find the mechanisms associa- 
ted with the rapid development of the three-dimensional waves. Raetz (1959) 
and Craik (1971) suggested the existence of a resonant triad consisting of a 
plane wave and a pair of symmetrical oblique waves. Their studies suggested 
that, if the wavenumbers are such that the phase velocities are equal, then the 
interactions might be strong. This is attributed to the phase-coupled energy 
exchange among the wave components and to the possible transfer of energy from 
the mean flow to the waves in the critical layer. Though the proposed model 
could not quantitatively explain the observed phenomena, it revealed the exis- 
tence of a particularly strong, nonlinear resonance mechanism for the selective 
growth of a pair of oblique waves and established in principle the subharmonic 
route to boundary- layer transition. Smith & Stewart (1987) considered the 
resonant interactions at wavenumbers close to the lower branch of the neutral 
stability curve by starting from the unsteady, triple-deck equations where the 
critical layer has moved into the viscous wall layer and is passive. The 
linear secondary stability analysis of Herbert (1983, 1988) explained many fea- 
tures of the observed phenomena in the parametric resonance stage. Other theo- 
retical and numerical work on the subject is reviewed by Craik (1986); Stuart 
(1986); Herbert (1988); Nayfeh (1988); Fasel (1990); Saric (1990); and Kleiser 
& Zang (1991). These studies have resulted in considerable effort to establish 
the mechanisms and methods for quantitative analysis in this area. 


1.2. Scope of Present Work 

The present work is concerned with developing a first-principles theory 
for use in the study of the fully interactive, near-resonance triad in 
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boundary- layer transition. The triad consists of a two-dimensional TS wave at 
the fundamental frequency u and of a pair of three-dimensional oblique waves 
at the subharmonic frequency «/2 . The oblique waves make equal and opposite 
ancles, ±0, to the flow direction, and their common streamwise wavenumber is 
nearly half that of the plane wave. Observation of boundary- layer transition 
typically reveals harmonic time-dependent instability waves, which exhibit spa 
tial downstream growth that is described in the early stages by linear stabil 
ity theory. As such, we will here consider initially linear instability waves 
whose continued downstream spatial growth results in nonlinearity (as was first 
done by Goldstein, Durbin, & Leib, 1987). Here, the fully nonlinear interac- 
tions of the waves have been studied via matched asymptotic expansion in the 
nearly common, critical layer at the transverse position where the instability 
wave's phase velocity equals that of the mean flow. 

Critical-layer nonlinearity (CLN) has been successfully adapted for the 
study of the nonlinear development of stability waves in various situations. 

Mas 1 owe provided an excellent review of CLN in shear flows: more recently, 
Huerre (1987) considered the two-dimensional shear layer in the viscous, 
critical-layer regime where growth, or nonequilibrium, effects are negligible. 
Critical-layer nonlinearity in the nonequilibrium regime, where ij *■ 8 
are dominant, has been considered by Goldstein & Leib (1988, 1989); Goldstein 
& Hultgren (1988); and Goldstein & Choi (1989), who showed that nonequi librium 
fi e growth) effects must be accounted for in order to match the nonlinear 
solution onto a suitable linear instability wave in the upstream region, where 
the wave amplitudes are small. They used their analysis to study various 
phenomena in free shear flows, such as nonlinear roll-up effects, the nonlinear 
evolution of a single oblique wave on compressible and incompressible shear 
layers, the nonlinear interaction of a pair of oblique waves in a free shear 
layer, and the instabilities in supersonic mixing layers. 


For boundary- layer flows, Goldstein, Durbin, & Leib (1987) used CLN to 
study the nonlinear roll-up of vorticitv in the presence of adverse pressure 
gradients, and Goldstein & Wundrow (1990) used CLN to study nonlinear second 
mode instability in high Mach-number boundary layers. The role of CLN in 
boundary- layer transition is evident from experimental observations, such as 
those of Corke & Mangano (1989), whose work demonstrates that nonl inear l ties 
are concentrated in the critical layer. The numerical results of secondary 
stability analysis (Herbert, 1988) also suggest that secondary instability on 
ginates from the redistribution of spanwise vorticity in the critical layer. 


Subharmonic resonance experiments in boundary layers (Corke & Mangano, 
1987, 1989; Kachanov & Levchenko. 1984; Saric & Thomas, 1983; and Sane et al . , 
1984) indicate that at least three stages of the nonlinear interactions follow 
the initial linear growth stage: super-exponential subharmonic growth, satura- 

tion, and decay. In the super-exponential stage, or, loosely termed, the para- 
metric resonance stage, the subharmonic oblique waves grow rapidly while the 
fundamental plane wave behaves almost in accordance with the linear theory. 

This stage is followed by the fully interactive regime, wherein the subharmonic 
exceeds the fundamental, saturates, and then decays. This regime, which may 
or may not be associated with a backreaction on the fundamental, ultimately 
results in the transition to the turbulent state. The recent experiments of 
Corke & Mangano (1989) and Corke (1989) have indicated the existence of several 
previously unencountered features of the fully interactive regime. Also, t e 
fully interactive regime has not been completely explored by direct numerical 
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simulations owing to the resolution requirements, which become prohibitive as 
the amplitudes increase downstream. For these reasons, this study focuses on 
the nonlinear interactions of the waves, starting from the linear stage and 
encompassing the resonance as well as the fully interactive saturation and 
decay stages. Emphasis is placed on capturing the fully interactive regime, 
which is done by allowing the oblique waves a sufficiently large amplitude to 
produce a backreaction on the plane wave. 

Because the aboye-cited experimental observations also indicate that sub- 
harmonic resonance occurs at low dimensionless frequencies, this study employs 
an asymptotic analysis that is valid not only in the low-frequency domain rele- 
vant to the experimentally observed phenomena but also at the technologically 
important lower- frequency domain, which may not be easily made subject to 
experimental testing or direct numerical simulations. 

Subharmonic resonance first occurs near, and continues downstream of, the 
upper-branch neutral stability curve (see the data of Kachanov & Levchenko, 
1984; Saric & Thomas, 1983; Saric, Kozlov, & Levchenko, 1984; Corke & Mangano, 
1987, 1989; and §3.4, part II of this study). The parametric resonance mech- 
anism is proportional to the plane wave's amplitude, which peaks at the upper 
branch; furthermore, as this study will demonstrate, the parametric resonance 
mechanism also is proportional to the cube of the Reynolds number, which 
increases indefinitely downstream. For all these experiments, these two fac- 
tors cause resonance to occur in the vicinity and downstream of the upper 
branch, with the subharmonic maximum growth rate occurring not at, but down- 
stream of, the upper branch. Therefore, the upper-branch scaling (Drazin & 
Reid, 1981) must be relevant to the subharmonic resonance phenomena. 

Moreover, at low frequencies, the upper-branch scaling applies not only in 
the vicinity of the upper branch, but also over most of the Reynolds -number 
range, breaking down only at low Reynolds numbers in the asymptotically small 
neighborhood of the lower branch (Goldstein & Durbin, 1986). Therefore, the 
upper-branch scaling is relevant to the entire Reynolds -number range associated 
with the subharmonic resonance phenomena. This scaling leads to a multideck 
stability structure characterized by a distinct, nonlinear, viscous critical 
layer separated from the viscous wall layer, and it cannot be captured by a 
triple deck structure (see Bodonyi & Smith, 1981; and Drazin & Reid, 1981, and 
references therein). 

The above factors uniquely fix the asymptotic scaling. As such, in this 
study a low-frequency, high-Reynolds-number scaling is adopted that is appli- 
cable to all the Reynolds-number range relevant to the observed subharmonic 
resonance. This scaling brings into play the critical layer and distinguishes 
it from the viscous wall layer. 

Section 2 of this study addresses the linear solution in the main part of 
the boundary layer. Slowly varying amplitude functions are introduced in the 
linear solution and are ultimately determined by considering the nonlinear vis- 
cous flow in the critical layer for each mode. The high-Reynolds-number linear 
solutions are obtained by matching the appropriate results in the different 
transverse zones. This leads to a linear relationship between the growth rate 
and the velocity jump across the critical layer. The waves then interact non- 
linearly, with the nonlinearity first coming into effect in the nearly common 
critical layer. 
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Section 3 addresses the nonlinear, viscous critical-layer flow. The non- 
linear viscous critical layer is required to approach a conventional linear, 
viscous critical layer far upstream. Within the order of approximation of the 
analysis, nonlinear effects are confined to the critical layer, thus providing 
the boundary conditions required to solve the viscous nonlinear equations that 
govern the critical-layer flow. The scales of amplitudes are determined for 
the fully interactive case, which implies that the oblique waves can become 
large enough to cause a backreaction on the plane wave. The self-interaction 
of the oblique waves occurs at the same level and is found to be of major 
importance. 

Section 4 addresses the nonlinear amplitude equations for the triad waves. 
The nonlinear solution presented here matches onto the upstream linear solution 
and onto the linear solution outside the critical layer. The amplitude equa- 
tions are obtained analytically by equating the velocity jumps across the crit 
ical layer that are calculated from the linear solution to those that are 
calculated from the nonlinear solution. 

Section 5 examines the nonlinear mechanisms independently of the flow 
divergence effects, which is feasible because of the comparatively short 
streamwise distance wherein nonlinearity occurs. Here, the study demonstrates 
that the saturation of the oblique waves is caused by their self- interact ion, 
even in the absence of the plane wave. This self-interaction, which leads to 
the saturation and decay of the oblique waves, previously has not been identi- 
fied in the subharmonic route to boundary- layer transition. The parametric 
resonance of the subharmonic caused by the plane wave, and the backreaction 
that occurs on the latter are also examined in §5. This is followed by the 
presentation of discussions and conclusions in §6. 

Part II of this study will examine the influence exerted over the develop- 
ment of instability waves by flow-divergence effects. This will be achieved by 
following Goldstein & Leib (1988) in forming a composite solution that accounts 
for both the nonlinear effects and for those caused by the boundary- layer 
growth. The results will be compared with experimental data as well as with 
numerical simulations, and excellent agreement will be demonstrated. 

Part II of this study also will examine other waves generated by the non- 
linear interactions, waves that were not present in the initial (upstream) 
linear stage and which are fully accounted for in the present theory.. These 
waves are of significant magnitude and play a crucial role in determining the 
triad's amplitudes as well as in interpreting several previously unexplained 
experimental observations. 

In addition, part II will examine the broad range of unstable spanwise 
wavenumbers . 


2. SCALING AND LINEAR SOLUTION 

The flow under consideration here is that of an incompressible laminar 
boundary layer. Its mean boundary- layer velocity is given by the Blasius 
velocity l T g, where 
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and X denotes the scaled Blasius skin friction « 0.332. The upstream flow 
starts as a near- resonance triad of spatially growing instability waves: a 

two-dimensional mode of normalized frequency u and wavenumber a, and a pair 
of subharmonic oblique waves of frequency «/ 2, streamwise wavenumber nearly 
equal to a/2, and spanwise wavenumber s#. The development of the waves up- 
stream of the nonlinear region (fig. 1(a)) follows the weakly nonparallel flow 
linear theory. All velocities are normalized by thg upstream velocity U^, 
lengths by 8*, time by S'/U^, and pressure by pl£, where p is the fluid 

density and 6* is the boundary- layer thickness defined as 8* --^/vx/U^. 

Here, v is the kinematic viscosity. , 

The high-Reynolds number scales as 

R.V°R- (22) 

where R is the Reynolds number based on the local-boundary layer thickness, 

8*, and R is an order-one scaled Reynolds number and a real quantity for 
time-periodic spatially growing disturbance. The small- frequency parameter a 
can be related to the normalized frequency F* by the approximate relation in 
Goldstein & Durbin (1988) 



The normalized complex wavenumber a is small, and its imaginary part is 
smaller than its real part (Goldstein, Durbin, & Leib, 1987). Consequently, 
each of the three modes has a well-defined critical layer at neaijly the same 
transverse position, y c , where the real part of their nearly common phase velo- 
city, c, is equal to the streamwise velocity, U. 


Outside the critical layer, the unsteady flow is governed by the linear 
dynamics, as first pointed out by Haberman (1972), and is treated as being 
locally parallel in the streamwise length scale over which the nonlinear 
effects take place. The velocity field is given by 


u 


U 0 (y) + &^eA Q (x 1 ) g— (y,x 1 ,o)e lX + 6A(x 1 )[u + (y,x 1 ,a)e lZ 

+ U_(y,x 1 ,a)e _lZ ]e 

v E Re i [eaA o <I> o e * X + SyA(e* Z + e* Z )$e* X ^j 


iX/2 


(2.4a) 

(2.4b) 


and 


w 
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E y 


(2.4c) 
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where 


(2.5a) 


_ 4 

Xj ; o X 


X = a(x - act) 
Z = [3z 


(2.5b) 

(2.5c) 


Y = 





(2.5d) 


The amplitudes A 0 and A 
respectively. The symbols 


are for the two-dimensional and the oblique waves, 
c and 6 are the measures of the amplitudes of 
the^twc^dimens i ona 1 "and” the oblique waves, respectively. The amplitude func- 
tions and their scales will ultimately be determined from the nonlinear ana y 
sis in the critical layer. Because the initial development of the instability 
waves is linear, we can take A(x x ) to be initially a real Quantity, but we 
allow the two-dimensional amplitude A 0 (xx) to be complex. Farther down 
stream, because of the nonlinear effects, both A(x a ) and , A o( x l) ^® come . c °J 
plex. In (2.4), U s and W* are related to the eigenfunction 4> through the 

relations 
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Y - H - 
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y ± 3y 


(2.6a) 
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(2.6b) 


1 JT-r 

where 4> n . * satisfy the linear Orr-Sommerfeld and Squire equations, respec- 
tively, with appropriate boundary conditions. The wavenumbers and phase velo- 
cities scale as (see Goldstein, Durbin, & Leib, 1987) 


a » aa + a 


4 A o 

'• 

(2.7a) 

a^ a A’ 

2i - A ’ 

(2.7b) 
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oc 

(2.7c) 

3 .. ’ 

a A 

ia A o 


ac 

(2 . 7d ) 

, a 3 A’ ’ 

2 - r 

ia 
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and 



a 


(2 . 7e ) 


where an overbar (~) denotes an order-one real constant. Here, (') denotes 
differentiation with respect to the flow variable xj. The scaling, given by 
equations ( 2 . 1 ) and (2.7), is consistent with Reid’s (1965) equation (3.128). 


At the high-Reynolds-number limit, the disturbance solution is a multi- 
zoned structure (fig. 1 (b)), wherein the critical layer is distinct from the 
wall layer (see Bodonyi & Smith, 1981; Drazin & Reid, 1981; Graebel, 1966; 
Eagles, 1969; Fraenkel, 1969; and DeVil lers, 1975). Zone I (Stokes layer) is 
a viscous wall layer and is the zone closest to the wall. Zone II (Tollmien 
region) is an inviscid rotational zone of adjustment within which the critical 
layer III is induced. Zone IV is above zone II and is an inviscid rotational 
region comprising most of the boundary layer. Above zone IV is the quasisteady 
zone V, in which the flow properties are of an inviscid irrotational type, 

U - 1 , V - 0. Zones IV and V can be combined into one zone (Heisenberg 
region) . 


The linear solutions in the above zones are presented in appendix A. 
Matching the inner solutions in zones II and IV produces relations (A. 4. 8 ) to 
(A. 4 . 11 ), which below are written for the lowest order of approximation as 

Y - a , (2.8a) 

X 3 -= 1 , (2.8b) 
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(*) 2 


cos © + 


'o (O 2 


A’ 


cos 0/ A 


4X 


3 (*♦) «■ ~T7= 

(c) ^/ 2 cRy 


8 \ 


U4> ) + 


-2 
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♦ iJUlf). 

0/ ini t i a 1 


Here, y c is the location of the critical layer, which scales as 


y - oY , 
J c c 


(2.9a) 


(2.9b) 


( 2 . 10 ) 


where Yc is an order-one real constant and and A$o are the velocity 
jumps across the critical layer for the oblique and two-dimensional modes, 
respectively. The obliqueness angle, 0 , is given by 


0 


sin 


-1 


L. 

Y 


( 2 . 11 a) 


8 



Upon substituting relations (2.5d) and (2.8a) into relation (2.11a), we 

0-60° . 

The velocity field in zone II, which contains the critical layer, 
according to the linear solution (appendix A) as 


u - oXY - a 4 f- (jr Y 2 ) + fie { [x + aa Q + a' 


‘( f ’ + ^cVo)] eA o eiX + 2 


+ c(X ± oa 1 tan 0sine + a 3 | f’ cos 0 + tan ^ ^XF(o) + V 


XY - c 


[(Y - Y e )ln|Y - Y c | * Y c In yj + XY 2 


1 V *| Y 

2 K C 4 ! \ C 2 


X Y 


XY - C 2 


4 - ~ sin 2 0 
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XY - c V XY - c cos^0/J 


+ + — r $ sin 2 ©! (cos Z)SAe* X ^ 2 /+..., 

cos e r xy - c 1 " 


*[ 


- 2 '" / /)„ [ : z-n , 2iy(cos Z)$Ae* X ^ 2 j + . . . , 


v = -o XY M iaeA e + 


w = -25(sin 0 sin Z) fie : v C * ~ iAe iX/2 + . . . , 

XY - C 

p - acX fie [cA Q e lX + 2 (cos 0 cos Z)6Ae lX/2 ] + . . . , 


where 


f - P, 


:( Y * v c ln ' V - Y J * Y c) * T y2 (? 


c Y 


and 


•“c ■ ■ 


X 2 . 2 


obtain 

(2.11b) 
is given 

r 

X + oa 
COS 0 


(2.12a) 

(2.12b) 
(2.12c) 
(2 . 12d) 

(2.12e) 
( 2 . 12 f ) 
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3. FULLY INTERACTIVE NONLINEAR SOLUTION IN CRITICAL LAYER 


The previous linear solution, (eq. (2.12)), becomes singular at the criti- 
cal layer. The appropriate transverse scaling coordinate in this region, which 
matches to the linear growth rate, is given by 



o 


(3.1) 


This scaling of the critical-layer thickness also is consistent with Reid's 
(1965) result that the critical layer has a thickness on the order of (aR) -1 ' 3 
Then, on substituting (3.1) into (2.12), the critical layer solution expands 
as follows 

o _ g . g i _o _2 —1 2 

u - c - o °Xn + o P c Y c n + o u c 2 n + So u_2 + So u_^ + Su^ + Sou^ + So Uj 


+ So 3 Ug + cu£ 0) + couj 0 ^ + eo 2 u^°^ + eo 3 Ug°^ + . . . , 


-eo _2 aXY fig iA e lX - Sa 2 yXY„ 2 cos Z fig iAe lx/2 + . . . , 


(3.2a) 

(3.2b) 


So 2 w 


-2 


So *w 


-1 


+ 5w r 


Sow 1 + 


So 2 w„ 


c 3 
So w„ 


(3.2c) 


p - eocX fig A Q e* X + SocX cos 9 2 cos Z fig Ae*^ 2 + . . . . (3. 2d) 


Here, u and u(°) denote the oblique and plane waves, respectively, and we 
have written 


v 


o 4 v 


and 



(3.3a) 


(3.3b) 


The full momentum equations can be expressed in terms of the scaled variables 
x^ , X , Z , and n as 


Du = -f aP Y + o 3 P v , o 8 P- , 0P 7 
~ \ a n 


and the continuity equation as 


au v + v- + (Jw 7 + 6 u - 0 , 

A f) Li 


(3.4a) 


(3.4b) 
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where 


u = (u,v,w) 


(3.4c) 


and 


. 3 -3 33 x 3 a 3 3 2 

D = a(u - c) «y + v + a u ~ + yi ~ — -2 ' 

8X 3n 9x i d£ R 3n 


(3.4d) 


It also is convenient to work with the equation for Z-component vorticity u, 
which can be written as 


afto - efa 2 u z w- + ouw z - o 6 av^ x - o 9 v z w x J 


(3.5) 


The nonlinear terms in the critical layer solution will balance the velo- 
city jump at the same order as in the linear solution when 


e - o 


10 


(3.6) 


The expansion (3.2) indicates that the first order of interaction occurs when 
S « e. This leads to a double exponential growth for the oblique waves, but 
the backreaction on the plane fundamental cancels out for this scaling level. 
The emphasis here, however, is on the fully interactive oblique and plane 
waves . 


The next order of interaction that produces a nonzero backreaction is 


7.5 


2.5 


(3.7) 


This also is the same order at which the oblique wave's self-interaction 
appears. Substituting (3.6) and (3.7) into (3.2), the expansions are now re- 
written in the general form as 


u - oc « g 3 Xn + g 5,5 Ug ^ + g^Ug + g^‘ 5 Ug 5 + 

5 . 5yr 6 W 6 • 5m /T^W X 

w * a W 5.5 + a W 6 + a *6.5 + a W 7 + " 

— 5 . 5,r 6 »t 6 . 5»t 7 y 

V ■ ° v 5.5 + ° v 6 + ° V 6 .5 + ° V 7 + ■ 

8.5 n 9 n 9.5 d ^ 10 D 

P “ ° P 8 .5 + ° P 9 + a P 9.5 + ° P 10 + 


and the vorticity, to, is given by 


-u = X + g 2,5 U, . - + g 3 U R - + g 3,5 U R - + g 4 U - + 
o .5 ,n 6 ,n 6 .5 ,n • >n 


(3.8a) 
(3.8b) 
(3.8c) 
( 3 . 8d ) 

(3.8e) 
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The mean flow components at the order a® and a® are given by 



“6 ‘ 

(3.9a) 


-2 

VI*U • 

(3.9b) 

The linear solution indicates that 



V c c - -c Re 2iy COS z Ae ix/2 , 
0 .o 

(3.10a) 


' : Y /2 

P„ _ - Xc Re 2 cos 0 cos z Ae , 

o . 5 

(3.10b) 


v° - -c Re iaA Q e iX » 

(3.10c) 


p u “ A o elX ' 

(3 .10d) 

Substituting 
each term in 

(3.8) into (3.4), we obtain the nonlinear viscous 
the expansion 

equations for 


L 3 U a “ " Yj L k u (3+a-k) " ^(3+1) ,x " p a,x ’ 

k=4 

(3.11a) 


L 3 W fi. " “ Yj L k W ( 3+2.— k) " ^ P (3+H),Z ’ 
k=4 

(3.11b) 


!!l (- 8 V g 3W * + 5 

3 - ’ -{* 3X 13 3Z 3x, J ’ 

(3.11c) 

where 

1 -It 3 1 3 2 

3 ' 3 811 ' R 3n 2 

/ 

(3. lid) 


L 4 ‘ c 3^ ’ 

(3. lie) 


L k = aU k 8X + V k ^ + ^ W k 8Z + U k-3 3x 1 ’ 

(3. Ilf) 
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u k- v k 


w k -0 


12 13 

2 ’ 2 ’ * ‘ ’ ’ 


U 3 ' * , 

(3 . llg) 

7 8 9 . 10 

for K - J. 2’ 2 • a ” d T ' 

(3 .llh) 


4. AMPLITUDE EQUATIONS 

In order for the nonlinear and linear solutions to match, the amplitude 
equations must be derived by equating , the velocity jump across the «itical 
layer obtained analytically from the nonlinear solution to that obtained from 
the linear one. The solutions at different levels in the expansion (3-? are 
needed to arrive at this result. The procedure is outlined in this section, 
and the details of the intermediate steps appear in appendixes B, C, and u. 

The interaction mechanisms contributing to the amplitude equation a PP®” in 
appendix E. The amplitudes for the triad waves appear here. °^ r "° n ^® arly 
generated waves, which determine the velocity jump calculated from 
linear solution and thus contribute to the amplitude equations of the triad, 
appear in appendixes B and D. These nonlinearly generated waves play a key 
role in interpreting several experimental observations, as will be discussed in 

part II of this study. 


4.1. 


Oblique Waves' Amplitude Equation 


As earlier noted, the amplitude equation for the oblique wave is obtained 
here by equating the velocity jump from the nonlinear solu ion with hatrom 
the linear solution (2.9a). Examination of the expansion (3.8) and the line 
solution (2.12) indicates that the velocity jump occurs at the U 1Q 5 leve . 
The equation for T,+ -*** rtKtain<lH frnm ,3 - 11 ‘ as 


(T „ - can be obtained from (3.11) as 
10 .5 ,n 


L 3 U l0.5,n " ^ W io.5 Yj [ L k u (i3.5-k),n + sin 6 L k,n W (i3 .5-k) , 

k=4 

- 3 COS e ( u k>z w (13 .5_k),n - U k,n W (l3.5-k),z)] 

where 

U + - U cos 0 + w sin 6 , 

— 3 3 

D + - cos 0 3 §2 " sin 6 a 8X ' 

The solution for U* , - can be written in the general form 

iu . o ,n 


(4. la) 


(4.1b) 

(4.1c) 
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u 


10. 5, n 


* E «-i> i[,n/2)Xt ” zl • 


n-0 ,m«l 


(4.2) 


Our purpose here is to match the velocity jump across the critical layer ob- 
tained from the nonlinear solution to that obtained from the linear solution. 

Since the latter is composed only of the oblique waves ~exp£i(|±Z^j , the solu- 
tion of - corresponding to n - 1, m - 1 is therefore needed. Analysis 

10 . 5)0 + 

of the Ul - equations indicates that U n - has no nm • 11 component for 
ft,n *>n 

ft - 5.5, 6.5, 7.5, 8.5, and 9.5. Thus, we can write 


L U +(11) - - 

^10. 5, n 


U c c -D + W_ + U Q -D + W, _ 
5.5,n a 8,n 5.5 


■ aU 5.5 U 8,nX " aU 8 U 5.5,nX ' K W 5.5 U 8,z)- 


- < W 8 U 5.5,z)- - V 5.5 U 


— - V„U* 


5.5 8,nn 8 5.5,nn 


(4.3) 


To obtain from the above equation, the solutions at the 5.5 and 8 levels 

must first be obtained. Using (3.11), the solution at the 5.5 level can be 
written as 

U 5 5 - 2 tan 0 cos Z Re ^iQ^n.x^e^^j (4.4a) 

W 5 5 - 2 sin Z Re Kn.xje^ 72 ] , (4.4b) 

where Q satisfies 


LjQ - yXcA sin 0 cos 0 , 


and using Fourier transform, we obtain as 

r° 

Q - A sin 0 


iKn Q hK /3 Av 
} e dK , 


(4.4c) 


(4 . 4d) 


where 



c 


(4.5a) 
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1 


(4.5b) 


f-XE 3 R 


x * x i - x 0 


(4.5C) 


and xg is the origin of the nonlinear region. 

The solutions for the 5.5- and 8-level components, provided in appendix B, 
are substituted into equation (4.3) to obtain Q{g ^5 ®fld it® Fourier transform 

(K). The velocity jump for the three-dimensional wave, J 3[) , is given by 

I 3 D - f «1 o.l d " - 9} o.l t |t - “> • (4 6) 


Details of the procedure are given in appendix C and the final result is 

3.5 


, i^A R p2 * • , u B 

X 2 ~ 4 0 V* 


2 * 
A^A , 


M = 1.055 - 17.8 i . 

The jump condition as given by equation (2.9a) is now written as 

f a 1 \ dA X 2 A . R T 

' cos e ' E V5 R x 30 

Substituting (4.7) into (4.8), we obtain the amplitude equation 


( a 1 X dA 
(, cos 6 * cos - ?/ E ■ 


which for 8 - 60‘ 


V2 R 
reduces to 


f *1) A + 5 ' t iA ’ A o 


p4 *5 2 

iM - A^A* 


X/X 


(4.7a) 

(4.7b) 


(4.8) 


(4.9) 


dA 4 k + 3 „ t iA - A ,| i,t! A V 

dx 5 o 10 X o 5 






Linear Resonance 


AaA 




interaction 


(4.10) 


where 
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(4.11) 


V X 2 2_E? 
0 “ ijl R " 8 X 3 

is the linear growth rate. 


4 . 2 . The Two-Dimensional Wave's Amplitude Equation 

Comparing the expansion (3.8) to the linear solution (2.12) indicates that 
the velocity jimp for the two-dimensional wave occurs at U 13 . The equation 
for U 13 is written from equation (3.11) as 


S u k ■ - Yj 

k-4 


which reduces to 
L 3 U 13 ,n * ^13 - ( L 5 


[ L k U (i 6 -k),n + sin 9 L k,n W (i 6 -k) 

- ® cos e ( U k,z W (l 6 -k),n ' U k,n W (16-k),z)] • 

.5 U 10.5,n + L 8 U 8 ,n + L 10.5 U 5.5,n) 

- sin e ( L 5 .5^ w 10 .5 + L 8 ,n W 8 + L 10.5,n W 5.s) 


( 4 . 12 ) 


+ B cos 01 




5.5,n W 10.5,z + 


U 8 ,n W 8 ,z 


IT _W 

10. 5, n 5.5, z 


U 5.5,z W 10.5,S- U 8,z*8,ii- U 10. 6 .z , 5.5 


,n) ' 


(4.13) 


Terms representing the interactions between components 6.5 and 9.5, 7.5 and 

8.5, or 6 and 10 would have produced an interaction term proportional to A . 
However, they failed to produce a velocity jump across the critical layer and 
as such are not written explicitly in equation (4.13). (Goldstein & Lee, 1991, 
observed the same phenomenon for the case of a triad of waves in an adverse- 
pressure-gradient boundary layer.) 
in the form 

u i3,n ’ 

n =0 ,m=l 

The linear solution at this level indicates that only the plane wave is discon- 
tinuous. Therefore, to match the nonlinear velocity jump to the linear one, 

only the plane-wave component q[ 20 ^ is needed. The procedure is similar to 

that outlined in §4.1; however, here the 5.5- and 8 -, and 10.5-level terms are 


Equation (4.13) provides a general solution 
q< 5”> (n,x 1 )e iI(n/2 ' X+mzl . (4.14) 
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needed to calculate The jump for the two-dimensional wave, J 2 q * * s 

then given by 


J2D-J’«‘3 0) - q‘ 2 3 0) (K - 0) . 


Details of the procedure are supplied in appendix D, which gives 


where 


i ffi R . . u r35 A o A * A „ a . a 3 

20 ‘ " r 5? A ° M i R ^ * " 2 x 


Mj - 0.5848 , 


M 2 - -0.241 . 


The jump condition (2.9b) can be written as 


dA_ 


a:' 


dx * 2f/l R 0 X J2D 


A„ - I 5 Jon + 


(4.15) 


(4.16) 

(4.17a) 

(4.17b) 


(4.18) 


^initial 

Substituting (4.16) into (4.18), we obtain the plane-wave amplitude equation as 


dx 


t k Q + ik t )A 0 + iMj 

X ^ ' X 

Linear 


d4.5 r® 3 

5 A AA* - iM„ % A*A d 

v/x 0 2 x 2 

J \ J 


Mutual 


Backreaction 


(4.19) 


where kj is the initial detuning factor. Equations (4.10) and (4.19) repre- 
sent the amplitude equations for the oblique and plane modes, respectively. 

The interaction mechanisms contributing to these equations are discussed in 
appendix E. Equations (4.10) to (4.19) indicate that the oblique subharmonic 
waves at 0 - 60° and the plane fundamental wave both reach the corresponding 

linear upper-branch neutral stability curve at R - 0.1537, which is obtained 
by setting the linear growth rate k 0 in equation (4.11) equal to zero. 


5. NONLINEAR MECHANISM 

Wave development is governed in part by the nonlinear dynamics considered 
above and by flow-divergence effects. Here, we will examine nonlinear mecha- 
nisms that occur over relatively short streamwise length scales where nonpar- 
allel effects are unimportant. The flow-divergence effects will be considered 
in part II, wherein a composite expansion will be obtined. 
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5.1. Oblique Wave Development in the Absence of the Plane Wave 

Because of the significance of oblique waves to boundary- layer transition, 
some experiments have attempted to examine oblique-wave development in the 
absence of the plane wave (e.g., Robey, 1987 and Schneider, 1989). This study 
examines the solution of the oblique waves' amplitude equation (4.9) in the 
absence of the plane wave. Under such conditions, the restriction on the span- 
wise wavenumber dictated by equation (2.8a) is lifted, hence the obliqueness 
angle can take any value. According to the asymptotic linear theory, the phase 
speed increases with the obliqueness angle and is given by 


c 


rX( 2 cos 0) 


(5.1a) 


For a general obliqueness angle, and with no plane wave present, equation (4.9) 
takes the form 


( cos 9 * 55Te) r 


X n 


R 2 


_ R rJl cos 0 - j-~* 

L a/ 2 4 X d (2 COS 0). 


R 4 ' 5 2 

A - iM(0) - A^A* 


XaA 


(5.1b) 


This equation is solved analytically to obtain 


A = Aj e 


q k K x 1+iC, 
M ob r »|» 


which gives the modulus as 


q k ,x 

I A| - A. e b G 


(5.2) 


(5.3) 


and the phase angle y as 

Y « C, an G . (5.4) 

The initial real amplitude at x = 0 is Aj , and G is the nonlinear function 
given by 


1 + 


Vi' 


- Vi' 


2 q 


k x\ - ^ ‘ ® 
ob x \ 


(5.5) 


where 
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(5.6a) 


d 


1 



5 


k ob 


I 


and the linear growth rate is given by 


\ 2 a/cos 6© ir R 

(5.6b) 

yR q X 3 cos 0 


1 

(5.6c) 

" " cos 0 + 1/cos 0 * 

r Real(M) 

i(r Imaginary(M) 

(5.6d) 


The case of 0 - 60°, relevant to the exact resonance case to be presented in 
part II, is discussed below. The effect of the obliqueness angle, correspond- 
ing to various spanwise wavenumbers, is discussed in §5. 1(b). 


5.1(a). Oblique Waves Self-Interactions at 0 - 60° 

The self-interaction of the oblique waves at 0 - 60° with no plane wave 
present is examined here. Equation (5.2) shows that there are only two parame- 
ters controlling the oblique-wave development. The first is the oblique waves 

amplitude, Aj , and the second is the scaled Reynolds number R. 

Figure 2 shows the modulus of the oblique waves' amplitude for several 
initial values, A* . The location with respect to the upper branch is given by 

r/r _ - 0.80, where R up - 0.1537 is the scaled Reynolds number at the upper- 
branch neutral stability curve. The initial growth depicted in figure 2 is 
linear It is followed by the oblique waves' nonlinear self-interaction, a 
mechanism that reduces their growth rate. Next, the amplitude of the oblique 
waves reaches an asymptotic saturation level independent of its initial level. 
This saturation level can be obtained from equation (5.3) as 


I A | 


'\nA k 


0.5 


ob 


M.R 4 ' 5 

l 


(5.7) 


Because the waves' self-interaction in equation (5.1) is proportional to their 
amplitudes cubed, a higher initial amplitude is associated with a quickening 
of nonlinear damping effects, which causes the asymptotic saturation level of 
the amplitude in the nonlinear region to be independent of the initial level. 

In the case of linear growth, the amplitude of the oblique waves remains 
real and its phase angle remains zero as initially prescribed. However, as 
figure 3 illustrates, nonlinearity causes the amplitude to become complex and 
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the phase angle to change with x. Increasing the initial level of amplitude 
for the oblique waves causes stronger nonlinear effects and rapid changes in 
the phase angle. In the limit of x ■* 00 , the phase angle decreases linearly 
with x at a constant rate independent of Aj , 

- I C k . . (5.8) 

di 5 * ob 

X-KD 


As earlier stated, the other factor controlling the self-interaction is 
the scaled Reynolds number. Figure 4 shows the linear and nonlinear solutions 

for two values of R corresponding tq R/R up * 0*5 and 0.95. At the smaller 
value of R, the linear growth rate is high and causes large amplitudes. Since 
the nonlinearity is proportional to A 3 , it produces a dramatic reduction in 

the amplitude. If R/R up - 095, the linear growth rate is weak, resulting in 
low amplitudes and less significant nonlinear effects. The net outcome is that 
the saturation level of the oblique waves slowly decreases with the Reynolds 
number, in accordance with equation (5.7). In reality, the Reynolds number is 
not frozen, as assumed here, but increases in the streamwise direction owing to 
the boundary- layer growth. Thus, the decrease of the saturation level with the 
Reynolds number (fig. 4) indicates that the oblique waves’ peak is followed by 
decay because the Reynolds number increases in the downstream direction. The 
oblique waves' self- interact ion thus not only causes their saturation, but 
ultimately also their decay. 


5.1(b). Effect of Obliqueness Angle 

The development of the oblique waves amplitude at various obliqueness 
angles corresponding to various spanwise wavenumbers is shown in figure 5. The 

initial conditions are R/R up = 0.8 and Aj - 0.1. Figure 5 depicts a linear 
growth region, followed by a nonlinear saturation region. The linear growth 
rate decreases with the increase of the obliqueness angle, as can be predicted 
from Squire's transformation. For all obliqueness angles, the self-interaction 
of the oblique waves leads to their saturation and ultimately to their decay. 

Of particular interest are the small-angle oblique waves, which are almost 
as probable as is the plane wave in a real flow with natural three dimension- 
ality. Figure 5 demonstrates that, at small angles, the linear growth is 
almost the same as that of the plane wave. Nonlinear effects are absent for 
the exact plane wave, 0=0°. But a small obliqueness angle results in pro- 
nounced nonlinear effects. The oblique waves are characterized by a three- 
dimensional vorticity field; therefore, the nonlinear effects appear at a lower 
amplitude level for the oblique waves than they do for the corresponding two- 
dimensional waves. Thus, for the present amplitude level, self-interaction is 
evident among the small-angle oblique waves but absent for the exact plane 
wave, 0-0°. This phenomenon is further depicted in figure 6, which shows the 
oblique waves’ saturation level at x « 100 versus the obliqueness angle. An 
abrupt drop in the saturation level is apparent at small obliqueness angles: 

In the intermediate range 10° < 0 < 65°, the dependence on the angle is weak; 
at higher angles, a sharp drop occurs. In a study by Schneider (1989) heating 
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elements were used to force oblique waves into a boundary layer. The develop- 
ment of the oblique waves was then followed by measuring the shear stress. The 
saturation of the measured oblique waves versus the obliqueness angle was qual- 
itatively the same as obtained here. 


5.2. Presence of the Plane Wave 

We will now examine the simultaneous development of the plane and oblique 
waves while the Reynolds number is frozen. The development of the waves at 
0 - 60° is described by equations (4.10) and (4.19), which are subject to the 

following initial conditions: amplitude of the oblique wave, jA^; amplitude 

of the plane wave, |A Qi |; 6111(1 ini tial , phase of the plane wave, t oi . The ini- 
tial phase angle of the oblique wave is taken to be zero. The initial detuning 
factor, kj, is taken to be zero and «|»oi * 3ir/2 * 

Figure 7 shows the development of the oblique waves for the fully interac- 
tive case. The initial conditions are R/R U p - 0.8, AqJ - 1, and Aj - 0.0001. 
In figure 7(a), the oblique waves' amplitude is compared with the amplitude ot 
the linear-growth case and with that of the parametric resonance case (only the 
linear and resonance terms are kept in (4.10)). The initial growth of the amp- 
litude is identical to the linear one; however, the first stage of nonlinear 
interaction, caused by the parametric resonance mechanism (the second term in 
eq. (4.10)), soon comes into effect, causing exponent i a 1-of-exponent lal growth. 
As the amplitude increases, a second stage of interaction (the last term in 
eq. (4.10)) comes into effect wherein the self-interaction of the oblique waves 
reduces their growth rate and leads to their saturation. 

In the linear-growth regime, the phase of the oblique waves amplitude 
remains zero. In the resonance regime, if toi “ 3w/2, tlie phase still remains 
zero, as depicted in figurji 7(b). As the fully interactive stages occur, the 
phase decreases with x, indicating a nonlinear influence on the streamwise 
wavenumber . 

Figure 8 depicts the development of the plane wave under the same initial 
conditions as in figure 7. In the fully interactive stage, the amplitude of 
the plane wave is reduced compared to that of the linear growth (fig- 8(a)). 
However, the effect on the plane wave of the interactions is not as pronounced 
as in the case of the oblique waves. In the plane-wave equation (4 19), the 
coefficient of the A 0 AA* term is purely imaginary and, consequently, merely 
results in a nonlinear detuning of the streamwise wavenumber. Also, the back- 
reaction term is a quartic-tvpe , A 4 . Thus, the effect of nonlinearity on the 
plane wave is felt farther downstream, when the oblique waves reach a consider- 
able amplitude. As for the phase of the plane-wave amplitude, figure 8(b) 
shows that it remains equal to its initial value up to the second stage of 
interaction, when it decreases in the downstream direction. 

Figures 7 and 8 show that the development of the waves can be divided into 
distinct stages, according to the mechanisms involved. In the initial stage, 
the linear-growth mechanism governs. In the following, the parametric reso- 
nance mechanism governs. The plane wave still follows the linear theory but 
causes exponential-of-exponential growth of the oblique waves. As the oblique 
waves continue to grow, the fully interactive regime comes into effect. The 
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oblique waves react back on the fundamental, causing a reduction in its growth 
rate. The most important mechanism in the fully interactive regime is the 
oblique waves' self-interaction. 


6. DISCUSSION AND CONCLUSIONS 

The effects of critical-layer nonlinearity (CLN) on a near-resonance, spa- 
tially growing triad of instability waves comprised of a plane wave at funda- 
mental frequency and a pair of symmetrical oblique waves at the subharmonic 
frequency were analyzed for a Blasius boundary layer. . A low-frequency, high- 

Reynolds-number scaling in terms of the small parameter o ~ leads to a 

stability structure with a distinct critical layer, wherein nonlinearity is 
important. 

In the fully interactive case, the amplitude of the oblique waves may 
exceed that of the plane wave such that a nonzero backreaction on the latter 
occurs. The amplitude equations were determined first by obtaining the analyt- 
ical solutions for each term in the asymptotic expansion of the nonlinear vis- 
cous critical-layer flow. Then, to match the solutions, the velocity jump 
across the critical layer, calculated from the nonlinear (inner) solution 
within the critical layer, was equated to the velocity jump across the critical 
layer calculated from the linear (outer) solution. Eliminating the velocity 
jump yielded a direct, nonlinear relationship between the growth rates and the 
amplitudes. 

The effort involved is justified by a comparatively simple amplitude equa- 
tion that differs from previously obtained equations. For one, the resulting 
theory is fully nonlinear in that all the nonlinearly generated oscillatory and 
nonosci 1 latory components produced by the interactions are accounted for. The 
analysis indicated that when the amplitude of the oblique waves is smaller or 
comparable to that of the plane wave, a parametric resonance term appears in 
the amplitude equation for the subharmonsc oblique waves. But the backreaction 
on the plane fundamental wave cancels out, even when the amplitudes are of the 
same order. A new stage of interaction comes into effect when the amplitude of 
the oblique waves becomes greater than that of the plane wave divided by o 2, °; 
the growth rates become fully coupled with a nonzero backreaction term affect- 
ing the plane fundamental wave. This is also the same order of magnitude at 
which the oblique waves self-interaction comes into effect. 

The development of the oblique waves in the absence of the plane wave was 
first examined in §5.1, where an analytical solution (5.3) was presented for 
the amplitude equation. This type of interaction in free shear flow was first 
considered by Goldstein & Choi (1989). The solution indicates that the sat- 
uration and decay of the oblique waves occur as a result of their self- 
interaction, a mechanism that does not require the presence of a plane wave or 
the appearance of a staggered pattern and that has not previously been identi- 
fied in the subharmonic route to boundary- layer transition. 

The study also indicated that, because the oblique waves initially follow 
the linear theory, their amplitude at first remains real, as originally pre- 
scribed. However, as the nonlinear effects become increasingly important down- 
stream, the amplitude changes from real to complex, indicating a nonlinear 
modification of the streamwise wavenumber. Equilibrium occurs on the short 
nonlinear length scale when the linear-growth mechanism is balanced by the 
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self-interaction mechanism.. The saturation level of the oblique waves is inde 
pendent of their initial level but would slowly decrease on the long length 
scale on which the Reynolds number increases. 

The saturation level of oblique waves was found to decrease with an 
increase in the obliqueness angle. As in the experimental data of Schneider 
(1989), increasing the obliqueness angle from zero to a small value resulted 
in an abrupt change in the nonlinear breakdown mechanism. The critical-layer 
singularity of the oblique mode was found to be stronger than that of the 
plane mode, which caused the nonlinearity and saturation of the oblique waves 
to occur at a much smaller amplitude compared to that of the strictly two- 
dimensional wave (of zero obliqueness angle). The saturation and decay mech- 
anism of the oblique modes was therefore found to differ from that of the 
plane wave. This result is in agreement with that found by Goldstein & Choi 
(1989) for oblique waves in free shear flows where the nonlinear interaction 
leads to explosive growth. 

Section 5.2 presented the results for the nonlinear mechanisms in the 
presence of the plane wave. The nonlinear effects were found to appear in two 
stages. In the first, the presence of the plane wave caused exponent ial-of- 
exponent ial growth of the oblique waves as a result of the parametric reso- 
nance mechanism. The oblique waves did not react back on the plane wave, 
allowing the latter to develop according to the linear theory, even when the 
oblique waves attained amplitudes on the same order of magnitude as the ampli- 
tude of the plane wave. However, as the oblique waves continued to grow and 

exceeded a certain level compared to that of the plane wave, a second fully 
interactive stage occurred. The growth rates became fully coupled and the 
oblique waves saturated because of their self-interaction while reacting back 
on the plane wave. 

Part II of this study will present results for the composite solution, 
which accounts for both the flow-divergence and the nonlinear effects. 

Therein, excellent agreement with observations is demonstrated not only over 
the parametric resonance stage, but also over the fully interactive saturation 

and decay stages. Novel features of the phenomena will be discussed and new 

interpretations of the experimental data will be given. 
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APPENDIX A 


LINEAR SOLUTION 

We present here the details of the linear solution in the five-zoned 
structure. 


A . 1 . Solution in Viscous Sublayer (Zone I) 

For the viscous effects to come into play in this region, the scaling for 
y should be 

y- o 4 y , (A. 1.1) 

where y is order one. This scaling is consistent with the classical theory 

result that the thickness of the wall layer is of the order (<xRC) 1/2 (Reid, 
1965). Substituting equation (2.1) into the Orr-Sommerfeld equation, we ob- 
tain, up to the required order of approximation, 

[ A J A A A A Art! A 

D 4 - iaR(Xo y - c)D^j* - 0 , (A. 1.2) 

AS A, _ _ ** _ 

where $ stands for $ or $, o stands for a or y, c stands for c or 

AS O As ~ U 

c, and D is the derivative with respect to y. We expand as 

^ 2 3 

* - g 0 + og 1 + a g 2 + a g 3 + . . . (A. 1.3) 

A As 

and, with a and c expanded according to equation (2.7), we obtain 

(8 4 + iR2c o D 2 )g o - 0 . (A. 1.4) 

The solution of (A. 1.4) must decay as y - « and is subject to the conditions 

g o (0) - g Q (0) - 0 . (A. 1.5) 

The first-order solution is thus given by 

s 0 ■ r (w - 1 ♦ .-*) . U.1.6) 

m 2 » -iaRc o , (A. 1.7) 
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where b Q is an arbitrary constant to be determined by matching with the solu 
tion at zone II. A similar solution is obtained for the oblique waves. The 
corresponding velocity components are given by 

U - /?« [cA Q b 0 (l - e -my )e lX + 2SAb cos 6 cos z(l - e _My )e lX/2 + . . .] 

(A. 1.8) 


A 


-io 5 a ~ (my + e~ my - l|eA o e iX , 


. 5 - b 


,-My 


- 2ia* y g iMy + e- , " jr - 1 |8A cos Z e *' + 


'- 1 ) 


iX/2 


(A. 1.9) 


and 


w - fie [2ib sin Z sin 6^1 - Me ** y )$Ae* X/r2 +•••]* (A. 1.10) 


where 



(A. 1.11) 


A. 2. The Tollmien Region (zone II) 

As pointed out by Goldstein, Durbin, & Leib (1987), the solution for this 
region is obtained by introducing the scaled transverse coordinate 

Y- y (A. 2.1) 

o 

directly into the Orr-Somraerfeld equation and using equations (2.1) and (2.7). 
The solution is obtained in the form 

<t> - oj^Cj + oc 2 + o 2 c 3 ) + (c 2 + oc 4 + o 2 c 6 )y] + a4 ^ • (A. 2. 2) 

As will be shown, the matching requirements give cj - C2 - C3 - 0; therefore, 
the solution can be written as 

$ - o(X + oa)Y + o^F(Y,$) , (A. 2. 3) 

where a is an order-one constant that depends on a, and F satisfies the 
relation / 




- \ X 2 Y(Y c + Y) 


(A. 2. 4 ) 
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where 


P 


c 


x 2 y 2 

C 


(A. 2. 5) 


and c - XY C . 

Since equation (A. 2 . 4) is singular at Y - Y c , F can be discontinuous 
across Y c , and F* denotes the solution above or under this point. Integra- 
ting equation (A.2.4), we obtain 

F*(Y,*) - F(0) + Y 2 * Y C [(Y - Y c )(ln|Y - Y c | + i* 1 ) 

* Y c (ln Y c ♦ ii')]} - V S (Y C ♦ Jy) . (A.2.6) 

where the normal velocity is continuous across the critical layer. The con- 
stants of integration and 4 >* are, in general, a complex function of x*. 

The value of F as Y - 0 is F( 0 ), which is determined via matching with the 
solution in zone I. The corresponding velocity components are given by 

u - Re eA 0 [(X + oa) + o 3 F(Y,4>)] e lX + 26A cos z[(X + aa)cos 0 + 


x sin 9 tan 0 a(X + oa)Y + a cos 0 F f ( Y 5 4 > ) 


A | T * 

+ o* n sin 0 tan 0 

U - c 


FJY^jje* 


iX/2 


and 


w «&^ 2 i SA sin Z sin 0^(X + oa) - ^ 3 c o(X + oa)Y 


+ a 3 F{Y,+) - ° 4 F(V 


,*)]e iX/2 } . 


(A. 2 . 7) 


^-eA Q iao 2 ^(X + oa Q )Y + o 3 F(Y,4i)j e* X 

- 2SAiyo 2 [(X + oa)Y + o 3 F(Y,$)]cos. Z e lX/2 j , (A. 2. 


8 ) 


(A. 2. 9) 
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A. 3. Solution in the Heisenberg Region (Zones IV and V) 

According to Goldstein, Durbin, & Leib (1987), Miles' (1962) solution, 
which is uniformly valid for y - 0(1) and y » 1 in the limit as a,y "* 0, is 
given by 


where 


D$_ U' _ 1 

^ ^ ^ 9 

$ U - c (U - c) Q* 


+ 0(a 5 ) 


a(l - c) 


+ aflj + a 2 0 2 + • ■ 


(A. 3.1) 


(A. 3. 2) 


. 1 (U - c) (1 - c), 

° (1 - C) 2 1(1 - C) Z (U - C V 


(A. 3. 3) 


1 (i - c y 


^2 (U - c)"Q 0 dy , 


(A. 3. 4) 


~ 9 2 

(U - c) 2 ^2 + 0 o dy 


(1 - c) 


(A. 3. 5) 


By substituting equations (A. 3.1) to (A. 3. 4) into the classical inv|scid 
function" Lin, 1955, p. 37), which is defined as 


- cD4> 

U'i - (U - c)D$ 


(A. 3. 6) 


and by inserting U = Uq and equation (2.7) into the result, expanding for 
small o, and then using equation (2.1), we obtain (following Goldstein, 

Durbin, & Leib, 1987) 

W = W + (y;a,5,9) + a 3 ^(cos 0 + ^-§) j~ , (A. 3. 7) 


iu 3 2icX “o 

W - W (y;o,c ,a) + a ’ 

a o 
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where 

W f (y ;o,C ,a) - ^(1 - aE )- 2 ° cX 


(1 - oc) 4 


[i, * *«J 2 + ° 2 <?h ‘s^ °(y 3 )] 


2 ~- 

+ a acX 


J + acj, --iy + °(y 2 ) 
4 5 4X J 


+ a 2 c 


5 lny '^> + ° (y) 


The coefficients J 1 to J 5 in the above equation are given by 


(A. 3. 9) 


n® 




n 2 1 _ 1 
U B " „2 + „_.*2 


Ug J 


dy , 


(A. 3.10a) 


'2 E 


u: 


h + u b 

U B 


1 2 
+ - 


(Xy) 3 T (Xy) 2 


dy , (A. 3. 10b) 



»® 

2 8 6 3 

8 ... + 6 f 1 


0 

"U 4 U 3 U 2 (Xy) 4 

(Xy) 3 (Xy) 2 " 4X 3 y(y + l) 


dy , (A. 3. 10c) 


J 4 = 2 


U‘ 


J o 


< - dy ■ 


0 

30 , 


J = lim -s- + X£2_ 

5 y-0 V3c V 


c+0 

u*u, 


B 


(A.3 . lOd) 


(A.3 . lOe) 
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A. 4. Matching Solutions 

Matching the solution in zone I and the solution in zone II is easily done 
by matching the velocities. Considering the^inner solution for the two- 

dimensional wave ((A. 1.8) to (A. 1.10)), as y ♦ ®, the solution reduces to 

u - fie b Q e iX , (A. 4.1) 


v = fie [-iao 5 b o (y - i) + • ■ -] elX • (A. 4. 2) 

Considering the solution in zone II ((A. 2. 7) to (A. 2. 9)), as Y - 0, we have 

u - fie [x + oa + o 3 F'(0) + . . .]e* X » (A. 4. 3) 

v - fie [-io 5 i o F’(0) + . . .]e lX . (A. 4. 4) 

Matching the u-ve loci ties gives 

b - X , (A- 4. 5) 

o 


while matching the v-velocities gives the imaginary part of F(0) as 


F. (0) 



(A. 4. 6) 


Matching the veloci tiesof the oblique waves produces similar expressions, with 
c replacing c 0 and y replacing a. 

Matching the solution in zone II and the outer solution in zones IV and V 
is most easily done by using the inviscid function. The inviscid function for 
the solution in zone II is obtained by substituting equation (A. 2. 3) into 
equation (A. 3. 6) and re-expanding 


W 


u 

X 


3 

+ o 





- S x V< 2V c - Y > 


+ 0 (< 7 4 ) 


(A. 4. 7) 


for Y > Y c> where A<j> denotes either ($ - 4> + ) or ($ 0 ~ Matching equa 

tions (A. 3. 7) and (A. 3. 8) shows that 


A 1 

a + + 0(o 4 ) 

A o 


(A. 4. 8) 
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. c ->-2 - ocX 

X r-(l - oc) - 1 - 


2-2 


(1 - oc) 4 


+ 2ctcJ 2 + a" c^Jo ) + o“acX| 


J 3 ) + o 2 «cx(l 4 * «J 5 ) 


+ a 


~ir ln Y c - 


48.2 -2 aJ XY + 4 

X a o/ c 


0 (A. 4. 9) 


and 


( c 


cos 0 + 


1 \ A' 


- 2 - 
Y M, 


cos 9/ A 


-2 

C Y 

y r^) ♦ — 


- 2 c 
CX^ 


cXY c ^/^cRy 


(A. 4. 10) 


A ' a 2 y 
o c 


Y J*4> ) * 


A - 2 c o 

O 2c q X 2 ° 2C 


-2 / a : 

a + i M- 


o Y c X 


(A. 4. 11) 


initial 


Equations (A. 4. 8] and (A. 4. 9) are dispersion relations that determine a 
and y in terms of c. Since their coefficients all are real, they are con- 
sistent with our original assertion that a and c are real quantities. In 
fact, it follows from equations (2.5d) and (A. 4. 8) that 


|3 - Ap a + — o Z 9m J 2 + 0(a 4 ) . (A. 4. 12) 

2 V3 o 

This shows that Q and a satisfy the usual long-wavelength, small-growth- 
rate resonance condition to within the order_of the detuning. To a first- 
approximation, equation (A. 4. 9) shows that c 0 and a satisfy the usual long- 
wavelength, smal 1-growth-rate dispersion relation 



The imaginary part of the dispersion relations (i.e., the matching between 
zones II and III) produces the real parts of equations (A. 4. 10) and (A. 4. 11). 

For the initial linear growth, A<j> is real. The amplitude A is real, and Aq 
is initially complex, but its imaginary part remains equal to its initial 
value. Nonlinearity causes both amplitudes to become complex. Since the amp- 
litudes appear as order a 3 in the dispersion relations, the imaginary parts 
of A' /A or A 0 /A 0 are balanced by the imaginary part of the corresponding 

Therefore, equations (A. 4. 10) and (A. 4. 11) in their complex form are valid 
for both the linear and nonlinear growth regimes. They relate the slow growth 
rates of the instability waves A' /A and A 0 /A 0 to the phase jumps A$ and 
A4 >q across the critical layer. To determine these letter quantities, it is 
necessary to consider the flow in the critical layer. 

Equations (A. 2. 3), (A. 2. 6) to (A. 2. 9), (A. 4. 5), and (A. 4. 6) show that the 
flow field in zone II can be written as 
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+ 



cA o e 


iX 



XY - c 


+ c(X + oa l tgn Q S i n © + a 3 • cos © + 


tan 9 sin 0 


XY - c ^ 


fxF(O) + 




[(Y - Y c )ln|Y - Y c | + Y c In yJ ♦ XY 2 % \> c ~ \r 


Y c + ! Y 


x r 

XY - c 2 


)D 


XA' . 2 fl 
+ r — sin 0 
yiA 


XY 






W - c cos^Q/J 


+ 


cos 0 




(cos Z)6Ae 



+ . . ■ 


(A. 


v - -o 2 XY Re [iaeA Q e lX + 2iy(cos Z)6Ae lX/2 ] + . . . (A. 

w - -2S(sin 0 sin Z) Re cX - - iAe lX/2 + ■ . • (A 

Xl c 

! p - acX Re [ eA o e lX + 2(cos 0 cos Z)SAe iX/2 ] + • • ■ > ( A 

where 

f -H C (V*Y C ln|Y - T c | ♦ fj - £ * h) • <* 


oa 

0 


4.14) 

4.15) 

4.16) 

.4.17) 

.4.18) 
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APPENDIX B 


NONLINEAR SOLUTION AT THE s 8 LEVEL 
At the <? 8 level, equation (3.11) gives 

- I ^2 <° 2) - l * “W ' 


i - 1 5 L\i< 2 °) 


aXni - i ^ IUq “ w/ - 2|3i tan 0 Q* - cyA tan © 
R 3n / 


£_ 
R 3n 2 


iXni - i ^)u^ 22) - - f 5yA tan 0 


" S I? W 8° 2) * 2 * i( N* ' \ ~ c ^'% - k % ] ’ 

V<°?> - -2ij3W< 02) , 
v< 2 °) - -iiU* 20 * , 

o,n o 

u(°°) - 2Ug 02 ^ , 

II (2—2 ) |.(22) 

U 8 8 

sin 0 W^ 22 ^ - -cos 0 U< 22 > , 

U7 ( 2—2 ) «(22) 

Wg ""S’ 

y(l±l) ^ ^(20) _ H ( 00 ) ^ ^(l±l) ^ y ( 00 ) _ y (2±2 ) _ y(l±l) _ g 

8 8 8 8 8,n 8,n 8 , n 

The substitution (eq. ((4.5)) and Fourier transform, defined as 

CO 

F(k) - J e _ikn F(n)d n , 


F ^) - fcj 


e ikn F(k)dk , 


(B.l) 

(B .2 ) 

(B.3) 

(B.4) 

(B.5) 
(B .6) 
(B .7) 
(B.8) 
(B .9) 
(B.10) 
(B .11) 

(B.12) 

(B.13) 
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are used to obtain the solution for equations (C.l) to (C.6) as 


where 


and 



,,(02) _ 
U 8 

-dAA*(S 1 + S 2 ) + bj , 

(B.14) 

w(° 2 ) 

8 

- dAA*(<! 

: sin 2 0 Sg + - S 2 ) + bg , 

(8.15) 

,,(20) 

u 8 

- 2d — hA 2 ^2 sin 2 0 + S,.) + bg , 

a 

(8.16) 


„(22) 

u 8 

- d §- hA 2 S 5 + b 4 , 
a 

(8.17) 

v (20) 

8 

- 28 dcA 2 h^2 sin 2 0 Sg + S 7 ) + b g , 

(8.18) 

v (02) _ 
8 

2 Be dAA*( 

4 sin 2 e S 8 4 S, - S 10 ) 4 b 6 , 

(8.19) 



d- r % , 

aXch 

(B.20a) 



h - , 

f-*E 3 R 

(B .20b) 


S 1 ’ 

^ e 1 * 3 ' 3 dk , 

— ao 

(B .20c) 


S 2 = - 

% 

- ikn li. 3 /n 
e k e 1 * 73 dk , 

(B.20d) 


J -05 


S 


3 


[( h/ 3 )l 
e 

(k 3 .k 3 ) +i 

( k r k 2 )n] 

i 

( k l - k 2) 

CM 




(B.20e) 
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° [ikn+(h/6)k 3 -(h/6)(k 1+ k 2 ) 3 +(h/3)(k^+k^) , 

e L H^kj + 


[ik n +(h/6)(k 3 +k 3 )] 




k 2 - ^dkj dk 2 dk , 


(B.20f) 

(B.20g) 


rO 


1 [ikn+(h/6)k 3 +(h/3)(k 3 + k 3 )-(h/6)(k 1+ k 2 ) 3 ] 


k e 


nO 


^I e 


x H^kj + k 2 - k^dk 1 dk 2 dk , 

[ik n +(h/6)(k 3 +k 3 )] 


(B.20h) 


H(kj - k)dk x dk , (B.20i) 


,0 


[(h/3 )(k®+k|)+t(k 1 -k 2 )n] 

( k i - k 2 ) 3 


dkj dk 2 , 


(B.20 j ) 


S 


9 



hk 3 /3 

e 


dk , 


(B.20k) 





hk 3 /3 

e 


dk , 


(B.20&) 


where H is the step function, and bi to bg are constants, determined 
through matching to the corresponding linear solution outside the critical 
layer. 
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APPENDIX C 


VELOCITY JUMP FOR THREE-DIMENSIONAL MODE 
Upon using equation (4.2) in equation (4.3), we obtain 


2 I 1 - z - EAyi [ 2 'V cos 6 * C0S ®( U 8 !w * U 8°^) - Sin ® W 8°Sn ] 

R on ) 


- £A >‘( U 8 “n * <s’) * 2 “ si " 9 *o 


+ 23 


« 2) )- 


cos 0 - sin 0 


Vs 


(02) 


* ^( U 8 ,22) + U 8 20 ’ C0S ®) 

- 2 si " 9 ^ V f) ■ (C1) 

The solution Q at the 5.5 level is given by (4.4), whereas that at the o 8 
level is given in appendix B. Taking the Fourier transform of equation (C.l) 
gives 


where 


'C.l - * hk3/3 ] k e ' 1 * 1 ' 3 (?0 ♦ q ♦ ♦ <>. ' 


y « 2 sin 0 » 

; n - 2 iA cos © + A„Q_ , 


(C. 2 ) 


- X 
a 


Xc 2 onn 


r - i — (3A cos 0 uj 02) + A sin 0 W* 02 * - A*U* 2 °* cos ©) , 

s i - ..-2 \ 8 ,nn 8 ,nn o.nn / 


a Xc 


’ 2 ixc 2 


26. r 2 ( QD (° 2 )) cos 9 - Vs 02) sin 0 + Qjuj 2 ® 1 cos 0 ] , 

n 

(«; n v 8° - V'S 2 ) • 


2 i sin 0 


aXc 3 


(C.3a) 

(C.3b) 

(C.3c) 

(C.3d) 


The velocity jump J 3 Q is obtained by 
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(C. 4 ) 


'3D ' f «} ”1 *> - ■ °> - I 0 + 'l * h * I 3 • 


where 


J; 


-hk 3 /3 

e 1 dk r 


j “ 0 1 1 1 2 ) 3 


J 0 


and 


- H 


] n - 2iriA cos 0 + 2 - S1 °.- AI 


'0 " - X 

a 


X 5 2 00 


Jl 2aX 2 c 3 h 


— —2 -2 

- 3 (A*I + AI„)A + ^-(A*I - AI 2 )A + 8 §3 I 3 - 4 ^ A *1 

a ay ay 


- ^A-AI, 




irX c h a y \ (3 (J 


5 w i ' 


2 -2 

2 sin z 9 (3 

.2-3 -2 

X c a 


4 A'I n * iAA- |I U - jj AI 13 - I 14 


’o ' 


-hk 3 /3~ .. 

e Qnn * ’ 


hk?/3 r ~ i 
e 1 k 1 [Q(-k 1 )Jdk 1 , 


hk 3 /3 
e 1 k„ 


Q'(-k 1 )]dk l , 


(C. 5 ) 

(C. 6 a) 

(C. 6 b) 

(C. 6 c) 

I 

(C. 6 d) 

(C. 7 a) 

(C. 7b) 
(C.7C) 
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,0 


hk:/3~ 


ILK.. /Or,* -1 

e 1 [QQM-k^jdkj , 


(C.7d) 


hk 3 /6 „ 
e 1 k 2 


-hk 3 /6,* 9 

e 2 Q 2 dk 2 dk x , 


J -k. 


(C.7e) 


,0 9 n® 

hkJ/6 9 

e 1 k 2 


-hk 9 /6 « 

e 2 k 2 Q dk 2 dk 1 , 


-k„ 


(C.7f) 



p° n 

p°° r 

! 6“- 

hk 3 /3 
e 1 kj 

ik n 

e Q 

«■ 

V 

— OD 

—00 


,* ^ ik 9 n 

i- (A*Q + AQ*)e dk dk dn , 
K 2 


(C.7g) 



nO „ 

r ® p 


hk^/3 

e 

ik.n 
e Qn 



n 


—00 ^ 

— CD 


m e ik2 ” dk 2 dk, dn , 
k 2 


(C . 7h) 



,0 „ 

«® r 


hk 3 /3 

e 

ik n 
e Q— 



n 

\ 

—00 ^ 

_ao U 


(VQ - AQ D e ik 2” a, dk, dn . <C.7i) 

k 2 2 1 



r ° 3 
hk 3 /3 

r ® 

ik n 

(hk 3 /e)+ik 2 n 

* 

e 

e Q n 

e v 

V 

—00 

—00 

—00 


-hk 3 /6~ 9 

e 3 Q 2 dk 3 dk 2 dkj dn . 


(C . 7 j ) 



r° 3 

-hkr/3 

p® 

ik i n o 

fhk 3 /6)+ik 9 n 

*10 - 

e 

e V 

e K 

* 

-00 ^ 

si 

—00 

—00 


-hk 3 /6 ,* 
e 3 k 3 Q 


x dk 3 dk 2 dk 1 dn . 


(C.7k) 
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11 


' 

(k> 3 >/3. 

hkg/2 r 


e k 1 2j k 3 

6 k 3 S(k 3 - k l * k 2> 

V « 

—CO 

-00 ' 


-hk^/6~- 
e Q 2 


x dk^ dkg dkg d^ , 


12 


1 

(k 3 +k 8 )h/3 , 

'* hkg/2 


e V 1 2J k* 

®T^- 8(k 3 - k l * k 2> 

V 1 

-00 

-00 


-hk “/6 ~ 
e 4 k 4 Q 


x dk 4 dkg dkg dk^ , 


13 


nO 


k? (h/3)(k?+k 3 ) 

-3 e v 1 3y QQ*6(k 2 + kj + k 3 )dk 2 dk x dk 3 , 


—00 —00 


14 


,0 


6 (kg + kj + k 3 )[A*Q - AQ*]dkg dkj dkg , 


where Iq to I 14 are given by equation (C.7). Substituting Ig to 
into equation (C. 6 ), we obtain 


- n . A A* 

J Q - 2 iriA cos © - 4ir sin 0 


T„ - 


* 6 2 h 0 


T & A 2 A* 

J 1 ■ * - .2-3, 5/3 
a X c h 


[-(3 + 

r + t 

9 —2 2 

L V a J 

id 6 O 

ay J 


x 4irA 2 A* 3 6 

2 " .2-3, 5/3 -2-4 
X c n ay 


-ffiT? 2 T g + j*g T 1q 


_ _ .BtrAV. JL [_ 2T _ 4T 1 

3 .2-3, 5/3 -2-2 L 11 13j ' 

X c h ay 


The T's are the reduced form of the I integrals and are evaluated to 


(C.72) 

(C.7m) 

(C.7n) 

(C.70) 

J 14 

(C. 8 a) 

(C. 8 b) 

(C. 8 c) 

(C. 8 d) 

give 
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T 0 - 0.5, T 2 - 0.591, T 3 - 0.73, T ? - 1.03, T g - 0.0625, 1 1Q 


- 0.188, 


T 


11 


tri _. 
32’ 13 


tri 

4 


Using a -JxR, c -VrA, and cos 0 » 1/2, and substituting the 
values for the integrals (C.9) in (C.8) we obtain 


j - - - 1 * 

J 0 4 y2 4 O 


a 2 a*r 3,5 


'* ' 0 845 ^ 


i2**p3 -5 

I_ - (0.21 - 5.391) A * 

1 jsA 


I 


3 


-11.45i A 2 A* 


t± 

/\/\ 


Thus, we can write 


h * h + J 3 “ M 


a 2 a*r 3,5 


/A 


where 


M - 1.055 - 16. 8i . 

The jump for the three-dimensional mode is thus given by 


J 


3D 


-iir A - | ttR 2 A A* + M 
4X^ 


r 3 »5 o 
S A V 


A A 


T 7 . (C.9) 

numerical 

(C.lOa) 

(C.lOb) 

(C.lOc) 

(C.lOd) 


(C.lla) 


(C.llb) 


(C.12) 
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APPENDIX D 


VELOCITY JUMP FOR THE TWO-DIMENSIONAL MODE 

Substituting equation (4.14) into equation (4.13) and using the transfor- 
mation (eq. ((4.5), we obtain 


2ni - h — iz 
k 3n 


>!”' ■ t "" 


-8 + T 8-8 + T 5 .5-10.5^ ’ 


( 0 . 1 ) 


where T^-g is the interaction of the mean flow with the o8-level components, 
given by 


T^_g - 2iy c <x cos 0 cA q - aip c cn Ug ^ wa v - ~p c v .g 


.cn 2 ui 2 ?^ cos 0 - 2p„ cos 0 Vj 2 °) . (D.2a) 

The interaction among the 8-8 components is Tg_g 

T 8-8 ■ ^( 2 Si " 9 " 8 ° 2 ) u 8 ! n 2) - 2 C0S 9 U 8 2 ' 2 , U 8 “ ) - C0S 9 U 8 ° 2)U ^ 2) 


- I cos e u<°°> U < 29) ) * « 1^( E V' - vf) 

V 


- 2 COS 0 


u (2-2) y (02) 
s.nn 


c 2 


(D.2b) 


T 5 5 _ 10 5 is the interaction df the 5.5-level components with the 10.5-level 
components 


5 .5-10.5 


— Ay cos 0 uji 1 * 

1 10.5 ,nn 

,c . 




5 cos 9 si 


n 0^ 


4q(3 COS 6 u{”> - sin 9 wj”) ) 

L C 


^ sin 0 Q \\l}\ 
-2 Tin 10.5 

.c . 


(D.2c) 


Other terms appearing in equation (4.13) do not contribute to the (20) 
component . 

The equations for the U 8 . Wo, and V 8 components are given in appen- 
dix B. The solution at the 10.5 level is obtained by utilizing equations 
(3.11) and (4.5) to obtain 
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(D.3a) 


pi - h 


3L\j(H) . !!u_ 

10,5 iX5 


and 


ni - 


,11 


h iLVi”) . !l 

an 2 / 10 - 5 ixc ' 


nl0.5 


-iv6U +(11) 

1yCU 10.5 


where 


R u - -3fj + 5(5 f 2 - |5 tan 0 f 3 + aA o tan 0 f 4 + iy(3Af 5 - A*f 0 ) 


- i tan 0(f 7 - fg) + 2iv c yAc 2 f 10 + j a tan © y c c‘f n , 


l - 


- 2 , 


R . - f S l - 2" f 2 * 5f 3 + a o f 4 - r ( f 7 - f s> - f-**c E2f n 


♦ i/»f 9 * cos e — , 


and 


f l - Q* u a 20 ^ 


QUg 02 ^ 


f, - Qw 


( 02 ) 

8 


f 4 - v 


- C ■ 

f 6 ' C * U 8 22) ' 


f 7 ' V, V 


(02) 


(D.3b) 

(D.3c) 

(D.4a) 

(D.4b) 

(0.4c) 

(D.4d) 

(D.4e) 

(D.4f) 

(D.4g) 

(D.4h) 

(D - 4 i ) 
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(0.4 j ) 


f e ‘ V (20) • 


f - w< 02 > 

(D.4k) 

f io ■ n , 

(D.4i) 

f u ■ « 1 ’ 2 • 

(D.4m) 

f . U (22) 
*12 u 8,n 

(D.4n) 

The Fourier transform of equation (E.l) is given by 


<X> 

'(201 e 1 * 3 ' 6 i 
Q f' - ^ ' 1 «o * • 

ld 2aXc * 

J k 

(D.5) 

The velocity jump is given by 


00 

. f -hk 2 /6~ 

,2 »- 2 L e ’'o*/ 

J o 

(D.6a) 

where 


R Q “ T M-8 + T 8-8 + T 5 .5-10.5 ’ 

(D.6b) 

The jump J 2 D can be written as 


*2D " ^M-8 + ^8-8 + ^5 .5-10.5 ’ 

(D . 7 ) 

where 


r iLi \ m iil_ Aco ce 

J M-8 2 - X A o 4.2 V OS 

Y x 

(D.8) 


This jump results from the mean-flow interaction with the linear part of V^ 20 ^ 
component. The other terms in T^_g (eq. (5.5)) produce no velocity jump. The 
J g _g is given by 
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' 8-8 


where 


|-COS 0 — 


“ /ft 


A AA*R 45 2 1/3 I - — A 3 A* ^ R 5 I 82 

0 8 aA « 


(D.9a) 


81 


X dx e- x3/3 , 


0 


(D.9b) 


82 


- §[* 3 «H]- M 


XjH(x 2 


- x 1 )dx 1 dx 2 dx , (D.9c) 


and, substituting cos 9 - 1/2, we obtain 


_3 j 3 

J g 8 - 0.85674 0 /iA- S— I 81 - 0.808 R 5 I 82 . 


(D.9d) 


The first term in the above equation results from the term in Tg_ 8 involving 
U* 00 ) alone interacting with the linear part of V 30 . The second term in the 
above equation is the result of the term in Tg_ 8 involving U* Other 

terms in T„ Q do not contribute to the velocity jump. 

o-o 

The jump resulting from the 5.5-10.5 interactions is written as 

J. . in , * G + G, + G„ . (0 - 10) 

1 5 .5-10.5 a b c 

Jump G a corresponds to the first square-bracketed term in equation (D.2c), 
jump Gb corresonds to the second square-bracketed term in equation (D.2c), 
and jump G c corresponds to the third square-bracketed term in equation (5.7). 
The fourth square-bracketed term in equation (D.2c) produces no velocity jump. 
G is given by 


C a ■ tw iM 1 ! sin 9 T - ' 3 N T - 


(D.lla) 


°b 


6i 


/(J_\ sin 0 A*A J 
1 _3-4 7/3 

\aj \ C h 


ft 2 

2 sin 0 T, + 1 (3 cos 9 - sin 0)T, 


-2 

a 


bl 


b5 


(D.llb) 
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(D.llc) 


4irB 2 

-2 

a 


A*A 
X 3 c 4 h 7/3 


2B 2 | 2 0 T 

=2 s,n 8 T cl 
La 


+ sin 0(6 cos 0 - sin 0)T K 
2a 5C . 


A AA* sin 3 0 
4nfl o T 

2 3 5/3 4C ’ 
« X c h 


In the above equations, terms involving T. result from interaction with the 

q.yUO) constituent of the or |. Terms involving T 4 result from 

the interaction with the Q* constituent of the 10.5 components, and terms 

n (021 (02) 

involving T g result from the interaction with the or 

constituent of the 10.5 components. The jump resulting from the 5.5- to 10.5- 
level interactions is written as 


where 


and 


p3 . 5 p5 * * .3 

J c . in _ « N, A A* A + N_ K * * 

J 5.5-10.5 1 ^ O 2 X 


N - 1.48 T . - 2.228 T . - -0.2188 
1 a4 c4 


N 2 - -0.2338 T al ♦ 1.62 T a5 ♦ 2.1 T bl ♦ 1.026 ? b5 - 2.098 T cl 

- 2.93 T _ - -0.15 

CO 

The total jump for the two-dimensional wave is given by 


. MR 3 ' 5 + M- P- A*A 3 , 

8 x 2 0 1 ,/x 2 X 


where 


- 0.5848 


M 2 - -0.241 . 


(D.12) 


(D . 13a) 


(D.13b) 


(0.14) 


(D.15a) 

(D.15b) 
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APPENDIX £ 

INTERACTION MECHANISMS CONTRIBUTING TO AMPLITUDE EQUATIONS 

Equations (4.10) and (4.19) represent the development of the oblique and 
plane waves, respectively. The interaction mechanisms leading to these equa- 
tions will be discussed here. 


E.l. Oblique Waves 


Equation (4.10) indicates that there are three mechanisms governing the 
development of oblique waves. Its first term represents the usual linear 
growth; its second term represents the, fundamental -subharmonic interaction that 

results from the interaction of the v 20 component of the plane waves with the 

5.5 components of the oblique waves; its last term represents the oblique 
waves ’ self- interact ion . 


The oblique waves' self-interaction term results from the velocity jumps 
across the. critical layer, ] y J 2> and J 3> given in appendix C: L results 

from the interaction of the V^" 1 ^ component with U* and W* ; h 

results from the interaction of ‘the (5.5) (1_1) components with Wj° and from 

the interaction of the (5.5)* -11 ^ components with uj 20 *; and J 3 results from 

the interaction of the V Q terms produced by nonlinearity with U g g , that is, 

from the interaction of V< 02) with (5.5) (1_1) components and from the 

interactions of V< 20) with the (5.5) (_11) terms. (These interactions are 

O 

summarized in table I.) 


F.2 . Plane Wave 

The plane wave's nonlinear development is governed by the jump across 
the critical layer. For the sake of discussion, this jump was split into 
three separate mechanisms: Im-8> 1 8—8 » ^ J5.5-10.5- (These interactions 

are depicted in table II.) 

j 8 (eq. (D .8) represents the usual linear growth. This jump results 
from the linear part of the V 20 component interacting with the mean flow. 

The other 8 components interacting with the mean flow produce no jump. 

The second type of interaction, Jg_ g (eq. (D.9), results from the self- 
interaction of the 8 components. The 5.5 components first self-interact to 
produce several 8 components: , Vg 11111 ^, and Wg with nm - (02), (20), 

(22), (00), (2*2 ) , and (1±1). The Ug° component interacts with the linear 
part of the V 2 ^ component to produce the first term in J g_g ( et l- (5.30)). 
The U 2-2 interacts with the U* 02 * to produce the second term in Jg_g. 
Other 8 components produce no velocity jump. 
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The third nonlinear mechanism, J5. 5-10.5 (eq. (0-10)). results from the 
interaction between the 5.5 components and the 10.5 components. The U^g 5 
component results from the interaction between the 5.5 components and the 
8 components, which produces several terms in the U^g.S equation. Most of the 
terms in the Uig .5 equation (4.13) produce no velocity jump. Only three 
terms , 


«* U 8 20) ' [ V 5 1 5 ) ( U 8 02) “8 2 )]’ V i 20, «; ' 

interact with the 5.5 components to produce the velocity jump. The interac- 
tions of these three terms with V 5 5 produce G a (eq. (D.lla)); their inter- 
action with Q n produces G 5 (eq. (D.llb)); and their interaction with Q 
produces G c (eq. (D.llc)). 
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Table 1. - Oblique waves’ jump across 
the critical layer 
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Tao.e 2. - Jump across the critical layer 
(a) Resulting from mean flow 8-component interaction. 
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(b) Resulting from 8-component self-interactions 
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Table -. - Concluded. 


(c) Resulting from 5.5-10.5 interactions. 
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Figure 1.— Structure. 



Figure 2.— An^Ktude of the oblique waves under variable 
initial levels, in the absence of the plane waves. The 
oblique waves’ self-interaction leads to their saturation. 
Note that die saturation level of the oblique waves is 
independent of the intial amplitude. 



Figure 3.— The phase angle of the oblique waves' 
amplitude in the absence of the plane wave. 
Nonlinearity causes the initially real amplitude to 
become complex. 
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Figure 4.—' The effect of the scaled Reynolds number on 
the oblique waves' nonlinear development In the 
absence of the plane wave (Aj - 0.1). The saturation 
level of the oblique wave stawfy decreased with the 
Reynolds number. 
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